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ABSTRACT:In this paper we have investigated the stability nature of Hopf bifurcation in a two 
dimensional nonlinear differential equation.Interestingly, our considered model exhibits both supercritical 
and subcritical Hopf bifurcation for certain parameter values which marks the stability and instability of 
limit cycles created or destroyed in Hopf bifurcations respectively. We have used the Center manifold 
theorem and the technique of Normal forms in our investigation. 
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I. INTRODUCTION 

From historical point of view, search on the existence of periodic solution plays a fundamental role in 
the development of qualitative study of dynamical systems. There are many methods for locating the periodic 
solutions, for example the Poincare return map, Melnikov integral, bifurcation theory etc. In this paper we have 
usedCenter manifold theorem, Normal forms and Hopf bifurcation theorem to study the behaviour of the limit 
cycle. 

The term Hopf bifurcation (also sometimes called Poincare Andronov-Hopf bifurcation) refers to the 
local birth or death of a periodic solution (self-excited oscillation) from an equilibrium as a parameter crosses a 
critical value. In a differential equation a Hopf bifurcation typically occurs when a complex conjugate pair of 
eigenvalues of the linearized flow at a fixed point becomes purely imaginary.The uniqueness of such 
bifurcations lies in two aspects: unlike other common types of bifurcations (viz., pitchfork, saddle -node or 
transcritical) Hopf bifurcation cannot occur in one dimension. The minimum dimensionality has to be two. The 
other aspect is that Hopf bifurcation deals with birth or death of a limit cycle as and when it emanates from or 
shrinks onto a fixed point, the focus. Recently, Hopf bifurcation of some famous chaotic systems has been 
investigated and it is becoming one of the most active topics in the field of chaotic systems. 

Hopf bifurcation has played a pivotal role in the development of the theory of dynamical systems in 
different dimensions. Following Hopf soriginal work [1942], Hopf and generalized Hopf bifurcationshave 
beenextensively studied by many researchers [3, 10,15, 19,20, 21, 24,29, 32, 33, 37]. 

As is well known, Hopf bifurcation gives rise to limit cycles, which are typical oscillatory behaviors of 
many nonlinear systems in physical, social, economic, biological, and chemical fields. These oscillatory 
behaviorscan be beneficial in practical applications, such as in mixing, monitoring, and fault diagnosis in 
electromechanical systems[l 1,15,25]. Also, the properties of limit cycles are very useful in modern control 
engineering, such as auto tuning of PID controller [39] and process identification [36]. Early efforts in Hopf 
bifurcation control [38] focused on delaying the onset of this bifurcation[l] or stabilizing an existing bifurcation 
[40]. 

II. CENTRE MANIFOLD THEOREM AND ITS ROLE IN HOPF BIFURCATION 

The center manifold theorem in finite dimensions can be traced to the work of Pliss [28], Sositaisvili [31] and 
Kelley [22]. Additional valuable references are Guckenheimer and Holmes [15], Hassard, Kazarinoff, and Wan 
[16], Marsden and McCracken [24], Carr [7], Henry [18], Sijbrand [30], Wiggins [37] andPerko[27]. 
The center manifold theorem is a model reduction technique for determining the local asymptotic stability of an 
equilibrium of a dynamical system when its linear part is not hyperbolic. The overall system is asymptotically 
stable if and only if the Center manifold dynamics is asymptotically stable. This allows for a substantial 
reduction in the dimension of the system whose asymptotic stability must be checked. In fact, the center 
manifold theorem is used to reduce the system from N dimensions to 2 dimensions [17]. Moreover, the 
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Centermanifold and its dynamics need not be computed exactly; frequently, a low degree approximation is 
sufficientto determine its stability [4] . 
We consider a nonlinear system 

± = /(*), x e R n (1) 

Suppose that the system (1) has an equilibrium point x 0 at the parameter value \i = [i 0 such that 
/(*o) = 0. In order to study the behaviour of the system near x 0 we first linearise the system (l)atx 0 . The 
linearised system is 

x = Ax, x 0 G R n (2) 
Where A = Df(x 0 ) is the Jacobian matrix of / of order n x n. The system has invariant subspaces E s , 
E u , E c , corresponding to the span of the generalised eigenvectors, which in turn correspond to eigenvalues 
having negative real part, positive real part and zero real part respectively.The subspaces are so named because 
orbits starting in E s decay to zero as t tends to 00 , orbits starting in E u become unbounded as t tends to 00 and 
orbits starting in E c neither grow nor decay exponentially as t tends to 00 . Theoretically, it is already established 
that if we suppose E u = <p , then any orbit will rapidly decay to E c . Thus, if we are interested in long term 
behaviour (i.e. stability) we need only to investigate the system restricted to E c . 

The Hartman-GrobmanTheorem[21] says that in a neighbourhood of a hyperbolic critical point x 0 , the 
nonlinear system (1) is topologically conjugate to the linear system (2), in a neighbourhood of the origin. The 
Hartman-Grobman theorem therefore completely solves the problem of determining the stability and qualitative 
behaviour in a neighbourhood of a hyperbolic critical point. 

In case of non-hyperbolic critical point, the Hartman-Grobman Theorem is not applicable and its role is 
played by the centre manifold thorem. The center manifold theorem shows that the qualitative behaviour in a 
neighbourhood of a non-hyperbolic critical point x 0 of the nonlinear system (1) with x E R n is determined by its 
behaviour on the center manifold near x 0 . Since the center manifold is generally of smaller dimension than the 
system (1), this simplifies the problem of determining the stability and qualitative behaviour of the flow near a 
non-hyperbolic critical point of (1). 

III. CENTER MANIFOLD THEOREM FOR FLOWS 

The statement of the Center Manifold Theorem for Flows is as follows : 

Let / be a C r vector field on R n vanishing at the origin (/(0) = 0) and let A = D/(0). Divide the spectrum of 
A into three parts, a s , a c , o u with 

Re A 

Let the (generalised) eigenspaces of a s ,G c anda u be E s , E c andE u , respectively. Then there exist C r 
stable and unstable invariant manifolds W s and W u tangent to £ s and E u at 0 and a C r_1 center manifold W c 
tangent to E c at 0. The manifolds W s , W u , W c are all invariant for the flow of /. The stable and unstable 
manifolds are unique, but W c need not be [7,15]. 
The system (1) can be written in diagonal form 

x = A c x + / 1 (x,y,z) 

y = A s y + f 2 (x,y,z) (3) 
z = A u z + h{x,y,z) 
with f t (0, 0, 0) = 0, Dfr (0,0,0) = 0 

/ 2 (0,0,0) = 0, D/ 2 (0,0,0) = 0 

/ 3 (0,0,0) = 0, D/ 3 (0,0,0) = 0 

Where/1,/2,/3 are some C r , (r > 2) in some neighbourhood of the origin,i4 c ; i4 s and A u on the blocks 
are in the canonical form whose diagonals contain the eigenvalues with Re A = 0, Re A < 0 and Re A > 0, 
respectively, (x,y,z) E R c x R s x R u , c = dim£" c since the system (3) has c eigenvalues with zero real part, 
s = dimE s and u = DimE u , f l9 f 2 and f 3 vanish along with their first partial derivatives at the origin. 
If we assume that the unstable manifold is empty, then (3) becomes 

x = A c x + /i(x,y) 

y = A s y + f 2 (x, y) (4) 

where A (0,0) = 0, £>A(0,0) = 0, 

/ 2 (0,0) = 0 D/ 2 (0,0) = 0, 
i4 c is a c X c matrix having eigenvalues with zero real parts, A s is an s X s matrix having eigenvalues with 
negative real parts, and f t and f 2 are C r functions (r > 2.) 

The following theorems established by Carr [7] in fact forms the basis of our discussion on the role of Centre 
manifold theorem in Hopf bifurcation and so we have stated them below : 
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Theorem3.1: 

There exists a C r center manifold, 

W c = {(x,y)|y = h(x), \x\ < S, /i(0) = 0,D/i(0) = 0}, 
for 8 sufficiently small, for (4) such that the dynamics of (4) restricted to the center manifold is given by the c- 
dimensional vector field 

u = A c u + f x (u, h(u)) (5) 

Theorem3.2: 

(i) Suppose the zero solution of (5) is stable (asymptotically stable) (unstable) then the zero solution of 
(4) is also stable (asymptotically stable) (unstable). 

(ii) Suppose the zero solution of (5) is stable. Then if (x(t),y(t)) is a solution of (4) with (x(0),y(0)) 
sufficiently small, there is a solution u(t) of (5) such that as t -> 00 

x(t) =u(t) + 0(e~ yt ) 
y{t) = h(u(t)) + 0{e-Y t ) 

wherey > 0 is a constant. 

From the theorem 2.2 it is clear that the dynamics of (5) near u = 0 determine the dynamics of (4) near (x, y) = 
(0,0). 

To calculate h(x) substitute y = /i(x) in the second component of (3) and using the chain rule, we obtain 

Af(/i(x)) = Dh(x)[A c x + f x {x, /i(x))] - Ch(x) - g{x, h(xj) = 0 (6) 
with boundary conditions /i(0) = D/i(0) = 0. This differential equation for h cannot be sloved exactly in most 
cases, but its solution can be approximated arbitrarily closely as a Taylor series at x = 0. 

Theorem 3.3 : 

If a function 0(x), with 0(0) = D (f)(0) = 0, can be found such that Af(0(x)) = 0(|x| p ) for some 
p > 1 as |x| -> 0 then it follows that 

h(x) = 0(x) + 0(|x| p )as|x| ^ 0 
This theorem allows us to compute the center manifold to any desired degree of accuracyby solving (6) to the 
same degree of accuracy. 

In the discussion above we have assumed the unstable manifold is empty at the bifurcation point. If we include 
the unstable manifold then we must deal with the system (3). In this case (x,y, z) = (0,0,0) is unstable due to the 
existence of a u-dimensional unstable manifold. However, much of the center manifold theory still applies, in 
particular Theorem 2.1 concerning existence, with the center manifold being locally represented by 

W c (0) = {(x,y,z) G R c XR S x R u \y = h x {x),z = /i 2 (x), /i;(0) = 0,D^(0) = 0,i = 1,2} 
forx sufficiently small. The vector field restricted to the center manifold isgiven by 

ii = Au + f(x, /i x (u), h 2 (u)) 

IV. HOPF BIFURCATION AND NORMAL FORM 

One of the basictools in the study of dynamical behavior of a system governed by nonlinear differential 
equations near a bifurcation point is the theory of normal forms. Normal form theory has been widely used in 
the study of nonlinear vector fields in order to simplify the analysis of the original system [5, 8,9, 14,15,23, 
26,47]. 

Several efficient methodologies for computing normal forms have been developed in the past decade 
[2,12, 13, 34,41, 42,43, 44,45,46]. 

The method of normal forms can be traced back to the Ph.D thesis of Poincare. The books by Van der 
Meer [35] and Bryuno [6] give valuable historical background. 

The basic idea of normal form theory is to employ successive, near identity nonlinear transformations 
to eliminate the so called non-resonant nonlinear terms, and retaining the terms which can not be eliminated 
(called resonant terms) to form the normal form and which is sufficient for the study of qualitative behavior of 
the original system. 

The Local Center Manifold Theorem in the previous section showed us that, in a neighborhood of a 
nonhyperbolic critical point, determining the qualitative behavior of (1) could be reduced to the problem of 
determining the qualitative behavior of the nonlinear system 

x = A c x + F(x) (7) 
on the center manifold. Since the dimension of the center manifold is typically less than n, this simplifies the 
problem of determining the qualitative behavior of the system (1) near a nonhyperbolic critical point. However, 
analyzing this system still may be a difficult task. The normal form theory allows us to simplify the nonlinear 
part,F(x), of (7) in order to make this task as easy as possible. This is accomplished by making a nonlinear, 
analytic transformation (called near identity transformation) of coordinates of the form 
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x = y 4- /i(y),where/i(y) = 0(|y| 2 ) as |y| -> 0. (8) 
Suppose that Df [J (x 0 ) has two purely imaginary eigenvalues with the remaining n — 2 eigenvalues 
having nonzero real parts. We know that since the fixed point is not hyperbolic, the orbit structure of the 
linearised vector field near (x, fi) = (x 0 , [i 0 ) may reveal little (and, possibly, even incorrect) information 
concerning the nature of the orbit structure of the nonlinear vector field (1) near (x, fi) = (x 0 , jU 0 ).But by the 
center manifold theorem, we know that the orbit structure near (x, fi) = (x 0 , is determined by the vector 
field (1) restricted to the center manifold. 

On the center manifold the vector field (1) has the following form 



y) \ImXQi) ReXQi) )\yJ V/ 2 (x,y,z) 

Where/ 1 and f 2 are nonlinear in x and yand A (X), AQ*) are the eigenvalues of the vector field 
linearized about the fixed point at the origin. 
Now, if we denote the eigenvalue 

= a(ji) + ico(fi), (10) 
Then for our assumption of non hyperbolic nature we have 
a(0) = 0, a>(0)*0 (11) 

Now we transform the equation (9) to the normal form. The normal form is found to be 

x = a(»x - o)(ji)y + (a GO* - b(ji)y)(x 2 + y 2 ) + 0(|x| 5 , |y| 5 ) 

y = co(ji)x + a(ji)y + (bQi)x + a(/0y)(* 2 + y 2 ) + 0(|x| 5 , |y| 5 )(12) 
In polar coordinates the equation (12) can be written as 

f = a(/j.)r + a(/j.)r 3 + 0(r 5 ), 

0 = eoQi) + b(ji)r 2 + 0(r 4 ) (13) 
Since we are interested in the dynamics near \i = 0, therefore expanding in Taylors series the coefficients in 
(13) about \x = 0 , the equation (13) becomes 

r = a'(0)jur + a(0)r 3 + 0(/j. 2 r, fir 3 ,r 5 ), 

0 = a)(0) + a/CO)^ + /?(0)r 2 + 0(/i 2 ,Mr 2 ,r 4 ) (14) 
where " ' " denotes differentiation with respect to 0 and we have used the fact that a(0) = 0. 
Our goal is to understand the dynamics of (14) for small rand \i. This is accomplished in two steps. In the first 
step, weneglect the higher order terms of (14) to get a "truncated" normal form and in the second step we study 
the dynamics exhibited by the truncated normal form as it is well known that the dynamics exhibited by the 
truncated normal form are qualitatively unchanged when one considers the influence of the previously neglected 
higher order terms [37]. 

Now, neglecting the higher order terms in (14) gives 
r = &\ir + ar 3 , 

0 = co + c\x + br 2 (15) 
where for ease of notation, we have 

a'(0) = d,a(0) = a, oj(0) = co, o/(0) = c, 6(0) = b 

For — 00 < ^- < 0 and \i sufficiently small, the solution of (15) is given by 



(r(t),0(t)) = 



lid 



a 



bd\ 

co + I c — —J \x 



t + 0 



0 



which is a periodic orbit for (15) and the periodic orbit is asymptotically stable for a < 0 and unstable for a > 0 
[15, 37]. 

In practice, a is straight forward to calculate. It can be calculated simply by keeping track of the 
coefficients carefully in the normal form transformation in terms of the original vector field. The expression of a 
for a two dimensional system of the form 



(IH 



0 -6)^Q + /7(*,y)" 



<o) 0 /Vyy \g(x,y) 

with /(0) = 0 = flf(0) = 0, is found to be [15, 16, 24] 
1 1 

^ (fxxx fxyy 9xxy 9yyy ) ~^ C^ty (-/xx ^/y ) <^xy (^9xx 9yy ) fxx9xx fyy 9yy ) 

In fact, all the above discussions gets converged in the famous Hopf Bifurcation Theorem which states as 
follows : 

Let the eigenvalues of the linearized system of x = (x), x E R n , E R about the equilibrium point be 

given by = a^u) ± i/3(ii). Suppose further that for a certain value of ^say \i = fi 0 , the following 

conditions are satisfied: 



www.ijmer.com 



Hopf-Bifurcation Ina Two Dimensional Nonlinear Differential Equation 



2 



= d * 0 



I. a(0) = 0,/?(0) = w^0 

(transversality condition: the eigenvalues cross the imaginary axis with nonzero speed) 

3. a ^ 0, where 

1 1 

^ \fxxx fxyy 9xxy 9yyy ) (/^y (/o: /yy ) «9xy (^9xx 9yy ) fxxQxx fyy 9yy ) 

where f rv denotes— — (x n ,y n ),etc. 
(genericity condition) 

then there is a unique three-dimensional center manifold passing through (x 0 ,jU 0 ) m R n x ^ and a smooth 
system of coordinates (preserving the planes \i = const.) for which the Taylor expansion of degree 3 on the 
center manifold is given by (12). If a 0, there is a surface of periodic solutions in the center manifold which 
has quadratic tangency with the eigenspace of X (Mo)> ^(Mo ) agreeing to second order with the paraboloid^ = 
— (a/d)(x 2 + y 2 ). If a < 0, then these periodic solutions are stable limit cycles, while if a > 0, the periodic 
solutions are repelling (unstable limit cycles). 

V. DISCUSSION ON THE BASIS OF OUR CONSIDERED MODEL 

We discuss the following modeltaken from [15]. 

x = y 

y = ii 2 y + fi 3 x 2 + n 4 xy (16) 
In our analysis, we fix fi 3 and /i 4 for simplicity. It is seen that, with suitable rescaling and letting (x, y) -> 
(— x, — y), for any \i 3 , 0 the possible cases can be reduced to two: \i 3 = 1 and /i 4 = ±1. 
Case 1: [i 3 = 1, fi 4 = — 1 (Super-critical Hopf bifurcation) 

In this case the system (2) becomes 

x = y 

y = \i x + \i 2 y + x 2 - xy (17) 

The fixed points of (17) are given by 

(x,y) = {±4~ih> 0) 

which exists for \i x < 0. We denote both the fixed points as P 1 = (x + , 0) = (+V — Mi> 0)and P 2 = (x_ ; 0) = 

(— V~ Mi, 0), where x ± = ±V — Mi- 

The Jacobian matrix / for the linearized system is 

J = ( 2 ° 1 

\2x - y \i 2 - x, 
The Jacobian matrix / at the fixed point P 1 = (V~ Mi, 0) is 

_ / 0 1 

~ UT-Mi M2 - 4~lh. 

The eigenvalues are given by 

u 2 - V=Mi ± V (M2 - aFMi) 2 + SaFMi 
A — 



1 _ 2 2 _ 2 

We observe that/^ > 0 and A 2 < 0 for \i x < 0 and for all . 
Thus the fixed point P 1 = (V~~Mi> 0) is a saddle point. 

The Jacobian matrix / at the fixed point P 2 = (— V - Jh, 0)is 

0 1 

/ = 



The eigenvalues are given by X = 



M2 + V^Mi. 

M 2 1 ±V 0*2 +V~M l) 2 l 



2 

^ _ ^2+V-Ml+V (m2+aFmY) 2 -8V-Mi ^ X ~ M2+V-M1-V (M2+aFmT) 2 -8V-Mi 
1 _ 2 2 _ 2 

The fixed point P 2 is stable for \i 2 + V - Mi < 0 i.e. ^2 < — V - Mi and unstable for^ 2 > V - Mi- 
Next, we verify the conditions for Hopf bifurcations at \i 2 = — V - Mi- 

(1) . For /i 2 = — V - Mi ^ Mi < 0» the eigenvalues on (x_, 0) are given by 

X 1>2 = ±iV 2 V-Mi(= ±iV _2x ) 

(2) . Also we have (/?eA(M 2 )) I =1^0 
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(3). Stability of the Hopf bifurcation: 

To study the stability of the Hopf bifurcation, we change coordinates twice, first to bring the point (— V - fa, 0) 
to the origin and then to put the vector field into standard form. 
Letting x — x — x_and y — y, we obtain 

j) = Gx_ o)G) + G 2 -xy) (18) 
The coefficient matrix of the linear part is A = ^ ^ — — = ^ 

The eigenvalues are X = ±i->J— 2x_(= ±ia)) where co = V - 2x_ 

Therefore the eigenvectors corresponding to the eigenvalue X are (+^^~2^~ 
Then, using the linear transformation 

where7 = ^ ^ — is the matrix of real and imaginary parts of the eigenvectors of theeigenvaluesA = 

±ij—2x_ , we obtain the system with linear part in standard form as: 

-V"2x_\ ni\ (f(u,v) 

where the nonlinear terms are 

■t 

f(u, v) — - 2x v 2 — uv, g(u, v) = 0 

and we have the following: 

1 2 

fu ^ >fuu 0 >fuuu 0 , 1 ,f uvv 0 ,/y yj —2x ^ ^ ^ yj —2x ' 

^it 0 ->9uu ^ ->9uuv ^ ' «0v ^ ->Qvv ^ ->Qwv ^ ">9uv ^ 

Substituting these in the relation 

1 1 

^ T~7 \fuuu fuvv 9uuv 9vvv\ * r . . \fuv (fuu fvv) 9uv (Suu 9vv) fuuduu fvv9vv\ 




lj uuu j uuu cj uuu c/uuuj 1 66(J 



d 2 f 



1 1 

We get a = that is a = ^=<0 for fa < 0 [as f YV denotes - 1 — (0,0)1 

The bifurcation is therefore supercritical and we have a family of stable periodic orbits. Thus, we can 
conclude that the system undergoes supercritical Hopf bifurcation for parameter values given by fa = — V - fa- 
Below, we have given support to our above mentioned theoretical result through numerical simulation. In the 
first column of Figure 1, figures a (7), a (77) and a (7/7 ) shows the phase orbitand the vector field and the figures 
b(I), b(II) and b(III) shows the corresponding time series plot for our considered model for different parameter 
values and initial points. In a(I) and b(I), the parametervalues are taken to be fa = — l,fa = —1.1 which is just 
before the occurrence of supercritical Hopf bifurcation at fa = — l,fa = ~ 1 (according to our theoretical 
result). Here, we have considered the initial point asx 0 = 0.4, y 0 = —0.01. The figure a(I) shows that the origin 
is stable and the orbit spirals to it. The same conclusion is supported by the time series plot shown in figureb(I). 
For the next figures in a(II) and b(II), we have considered the parameter values fa = — 1, fa = —0.9 which is 
just after the occurrence of supercritical Hopf bifurcation. We considered the initial point as x 0 = — 0.66,y 0 = 
—0.06, which is inside the limit cycle produced during the Hopf bifurcation.The figure a(II) shows that the 
origin is unstable and the orbit spiral away from it onto the smooth invariant closed curve encircling it (the limit 
cycle). The same conclusion is supported by the time series plot shown in b(II). In figure a(III) and b(III) we 
have kept the parameter values same with the earlier case but changed the initial point as a point outside the 
limit cycle. Here we have considered fa = — l,fa = ~ 0.9 with initial point x 0 = 0.6, y 0 = O.Figure a(III) and 
b(III) shows that the orbit spiral towards the invariant circle. So, our numerical simulation clearly supports our 
conclusion that that the system undergoes supercritical Hopf bifurcation for parameter values given by fa = 
-yFfa. 
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Figl.The Phase orbit and time series of Eq. (16).a(I) and b(I) before the Hopf bifurcation at the parameter 
value p 1 = —1,P2 = ~ 1-1 w ith initial point x 0 = 0.4, y 0 = —0.01 , the origin is stable and the orbit spirals to 
it.a(II) and b(II)after the Hopf bifurcation at the parameter value p 1 = —1,^2 = —0.9 with initial pointx 0 = 
— 0.66, y 0 = —0.06, which is inside the limit cycle, the origin is unstable and the orbit spiral away from it and 
onto a smooth invariant closed curve encircling it.a(III) and b(III)after the Hopf bifurcation at the parameter 
value p 1 = —1,P2 = — 0.9 with initial point x 0 = 0.6, y 0 = 0 which is outside the invariant circle(limit cycle), 
the orbit spiral towards the invariant circle. 
Case 2.fi ± = 1, fi 2 = 1 (Subcritical Hopf bifurcation) 
In this case the system (16) becomes 

x = y 

y = p 1 + p 2 y + x 2 + xy (19) 

The fixed points of (19) are given by 

(x,y) = (±7^,0) 

which exists for ^ < 0. We denote both the fixed as P 1 = 0) = (+V — Pi> 0)and P 2 = (x_, 0) = 

(— V - Mil 0), wherex+ = ±V — Mi- 

The Jacobian matrix J for the linearized system is 

0 1 

/ = 



2x + y Jp~ 2 + 



x 
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( o 1 

The Jacobian matrix/ at the fixed point Pi = (\f-ju, 0)is / = n / , i 

_ _ M2 + V^Mi ± V (M2 + V^Mi) 2 + SV^Mi 

2 ; 

1 _ 2 2 _ 2 

We observe that A 1 > 0 and A 2 < 0 for ^ < 0 and for all \i 2 . 

Thus the fixed point P 1 = (V - Mi; 0) is a saddle point. 

/ 0 1 

The Jacobian matrix / at the fixed point P ? = (— V~~m7, 0) is / = 0 , > 

The eigenvalues are given by A = 

2 

^ ^ ^ _ ^2-y-^i+V (^2-aFmi) 2 -8v-^i ^ a - (M2-aFmi) 2 -8V-^1 

1 _ 2 2 ~~ 2 

The fixed point P 2 is stable for \i 2 — V - Mi < Oi.e. ^ 2 < V - Mi an d unstable for ^ 2 > V - Mi- 
Next, we verify the conditions for Hopf bifurcations at ^u 2 = V - Mi 

(1) . For jU 2 = V - Mi ? Mi < 0, the eigenvalues on (x_, 0) are given by 

A 1>2 = ±iV 2 V-Mi(= ±iV _2x -) 

(2) . Also we have observed that T~ ^(^2)) =1^0 

dfe '^2=^20 

(3) . Stability of the Hopf bifurcation: 

To study the stability of the Hopf bifurcation, we change coordinates twice, first to bring the point (— V - Mi> 0) 
to the origin and then to put the vector field into standard form. 
Letting x — x — (— V~ Mi)and y = y, we obtain 

,y) ~ (-2^^ 0) (y) + (x 2 + xy t 
The coefficient matrix of the linear part is A = (^_2^jzrj]~ q) = Gx 0) 
The eigenvalues are X = ±i A /2V— Mi = ±i^/— 2x_(= ±io)) where = V — 2x_ 
Therefore the eigenvectors corresponding to the eigenvalue A are (_|_j t; 2x 
Then, using the linear transformation 

where7 = ^ ^ — is the matrix of real and imaginary parts of the eigenvectors of the eigenvalues 

= ±i-yj—2x_ , we obtain the system with linear part in standard form: 

0 -V"2x_\ , U \ (f(u,v) 
0 J™* \9(u,vl 




where the nonlinear terms are 



1 2 



f (u, v) = UV + , V 

g(u,v) = 0 

and we have the following: 

1 2 

fu V , f uu 0 ,f uuu 0 , f uv 1 , fuvv 0 > fv U ~\- ^j _2 x 2^ 0 y fvv yj —2x ' 

Qu 0 > Quu 0 > Quuv 0 > Qv 0 > Qvv 0 > Qvvv 0 > Qu9v 0 

Substituting these in the relation 

1 1 

^ [fuuu fuvv Quuv 9vvv\ [fuv (.fuu fvv) 9uv (.Quu Qvv) fuu Quu fvvQvv\ 

1 1 d 2 f 

We get a = that is a = — ^=># for < 0 [as f YV denotes (0,0)1 

& 16x_ 16^J=i^[ ^ i L Jx y dxdy K ' JX 

The bifurcation is therefore subcritical and we have a family of unstable periodic orbits. Thus, we can conclude 
that the system undergoes subcritical Hopf bifurcation for parameter values given by \i 2 = V - Mi- 
Below, we have given support to our above mentioned theoretical result through numerical simulation. In the 
first column of Figure 2, figures a(I), a(II) and a(III) shows the phase orbit and the vector field and the figures 
b(I), b(II) and b(III) shows the corresponding time series plot for our considered model for different parameter 
values and initial points. In a(I) and b(I), the parameter values are taken to be = — \,\i 2 — 0.9 which is just 
before the occurrence of subcritical Hopf bifurcation at fa = —1, [i 2 = 1 (according to our theoretical result) 
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with initial point x 0 = —1.7, y 0 = 0 which is inside the limit cycle. The figure a(I) shows that the origin is 
stable and the orbit spirals to it. The same conclusion is supported by the time series plot shown in figure b(I). 
For the next figures in a(II) and b(II), we have considered the parameter values p 1 = —1, p 2 = —0.9 which is 
before the occurrence of subcritical Hopf bifurcation. We considered the initial point as x 0 = — 1.8,y 0 = 0, 
which is outside the limit cycle produced during the Hopf bifurcation. The figure a(II) shows that the origin is 
stable and the limit cycle is unstable. The orbit spiral away from the limit cycle. The same conclusion is 
supported by the time series plot shown in b(II). In figure a(III) and b(III) we have considered parameter values 
p 1 = — l,p 2 = 1-1 which is just after the subcritical Hopf bifurcation. The initial point was taken as x 0 = 
— 0.9, y 0 = 0. Figure a(III) and b(III) shows that the origin is unstable, the orbits spiral away from it and no 
limit cycle exists in the neighbourhood. So, our numerical simulation clearly supports our conclusion that that 
the system undergoes subcritical Hopf bifurcation for parameter values given by p 2 = V - Mi- 
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Fig2. The Phase orbit and time series of Eq. (4). a(I) and b(I) before the Hopf bifurcation at the parameter 
value p 1 = —1,1*2 = 0-9 with initial point x 0 = —1.7, y 0 = 0 , which is inside the limit cycle, the origin is 
stable and the limit cycle is unstable, the orbit spirals to the origin. a(II) and b(II) before the Hopf bifurcation 
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at the parameter value p 1 = —1, p 2 = — 0.9 with initial point x 0 = — 1.8, y 0 = 0, which is outside the limit 
cycle, the origin is stable and limit cycle is unstable and the orbit spiral away from it. (a(III) and b(III)) after 
the Hopf bifurcation at the parameter value p 1 = —1,^2 = 1-1 with initial point x 0 = — 0.9,y 0 = 0, the origin 
is unstable, the orbits spiral away from it and there is no limit cycle exists. 

VI CONCLUSION 

We have presented Hopf bifurcation of a nonlinear system, which admits both supercritical and 
subcritical cases depending upon different parameter values. The results were derived theoretically with the help 
of center manifold theory and normal form method. Numerical simulations were presented in support of the 
conclusions drawn theoretically. 
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